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On the lattice searching for the gluon condensate is difficult because a large perturbative con- 
tribution to the expectation value of the action has to be subtracted before looking for a small 
contribution from a possible gluon condensate. The perturbative calculation therefore has to be 
very precise. We use a modified version of stochastic perturbation theory to calculate a pertur- 
bative series in a boosted coupling, which converges more rapidly than the series with the usual 
lattice coupling, reducing the uncertainties in our results. We do not see any condensate of di- 
mension two, as suggested by some earlier lattice studies, but we do find a contribution from a 
dimension four condensate. The value of this condensate is approximately 0.04(1) GeV 4 , but 
with large uncertainties. 
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1. Introduction 

There is a long history of attempts to determine a gluon condensate on the lattice [[TJ, |2|] by 
looking at the expectation value of the gluonic action, which in the case of the Wilson action is 
simply the plaquette expectation value 

F = l(Tr(C/ D )) . (1.1) 

This is difficult, because the signal from any non-perturbative condensate will be much smaller 
than the background caused by perturbative gluon fluctuations, dominated by the contribution from 
gluons with momenta near the ultra-violet cut-off. The relationship between AP, the difference 
between the perturbative plaquette and the non-perturbative value measured in a Monte Carlo cal- 
culation, and the gluon condensate introduced in ^ is given by 



AP(g 2 ) = P per ,{g 2 )-PMC(8 2 ) = ^ 
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where a is the lattice spacing, j3(g) the j8 -function, and —bog 3 the first term in the expansion of 



In eq.(1.2) we make the conventional assumption that the condensate term has dimension four, 
so that AP will scale like a 4 . There has however been much interest in the possibility, suggested 
in that the dominant correction to perturbation theory scales like a 2 instead. 

In this work we concentrate on quenched SU(3) with the Wilson gauge action. Most of our 
results are based on 12 loop series found on a 12 4 lattice, though we will also present some ex- 
ploratory 16 loop perturbation theory results from 8 4 lattices. 

A long perturbation series can also be used to investigate other interesting issues. For example, 
how long does the series need to be before we see that we have an asymptotic series rather than a 
convergent one? 



2. The Perturbative Series 

In the introduction we have argued that a long perturbation series is needed to answer the ques- 
tions that interest us. In this section we discuss a number of methods of finding p n , the coefficients 
in the series for the plaquette, 

Ppert(g 2 ) = l+'EPng 2n ■ (2.1) 

2.1 Diagrammatic Lattice Perturbation Theory 

The first method that suggests itself is to use diagrammatic methods, similar to those used in 
the continuum. On the lattice the Feynman rules are more complicated than in the continuum, there 
are extra vertices, and the integrands in the loop integrals are much more complicated than in the 
continuum. Nevertheless the first three coefficients in the series have been found by conventional 
techniques [Q, [5j|. The advantage of this technique is that the coefficients can be found with very 
high accuracy, for example in [||] p2 is given to 9 significant figures, and to 6, which is much 
better than what can currently be achieved with stochastic methods. The disadvantage is that the 
difficulty of the calculation grows very rapidly with the number of loops, and it will probably be a 
long time before we gain many more terms in the series. 
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Figure 1: The first 16 coefficients in the perturbation series for the plaquette. Filled symbols show positive 
coefficients, empty symbols show negative coefficients. Errors are smaller than the symbols. 

2.2 Stochastic Perturbation Theory 

The solution to this problem is stochastic perturbation theory, a mechanised way to do lattice 
perturbation theory introduced in [f7|]. Each field of lattice gauge theory is replaced by a series in g, 
and the result of a simulation is a Monte Carlo estimate for each p n . The precision not as large as 
that found by doing integrals, but the effort needed grows as a polynomial of loop number, so we 
can go much farther in n. In [Q] the coefficients up to p% were found, [f|] extended this to the 10 
loop level, and in Fig.|l] we present exploratory results up to 16 loops, calculated on an 8 4 lattice. 

Unfortunately the series converges rather slowly, and simply truncating the sum at 16 loops 
does not give a good estimate of the complete sum. As can be seen from the graph, the dependence 
of the coefficients on n is quite smooth, so one possible approach is to fit a curve to the known 
coefficients and use this to estimate the p n for higher n. This was the approach we used in [S]. 
The new coefficients agree very well with the predictions made in [Q] when only 10 loops were 
available. Even at 16 loops we do not yet see signs of the series diverging, it still looks like a series 
with a finite radius of convergence. This suggests that the crossover region between the strong 
coupling and weak coupling regions is still having a strong influence on the series, as discussed 
in [Q]. Motivated by this behaviour in the perturbative series, the possibility that there is a third- 
order phase transition in this region was recently examined in jJTotl - 

2.3 Boosted Perturbation Theory 

It is well-known that perturbation theory in the lattice coupling g 2 converges slowly. The ex- 
planation is that the Aqcd parameter of the bare lattice coupling g 2 is surprisingly small, A/ a , = 
0.035A^ so g 2 is not a very good expansion parameter. The remedy is to use a "boosted" cou- 
pling which has a more reasonable Aqcd- The most popular choice is = g 2 /P which has 
A/a; = 03&Ajfs . The boosted coupling is larger than g 2 , but we hope that the coefficients will 
now decrease faster, and that we will find a reliable answer with fewer terms in the series. If we 
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Figure 2: The boosted perturbation series (black curve, square symbols) compared with the series in the 
lattice coupling (blue curve, round symbols). Filled symbols show positive coefficients, open symbols show 
negative. 

had the conventional series to very high accuracy, transforming to the new coupling would be easy. 
However if we only have limited accuracy in the coefficients the transformation is numerically dif- 
ficult as we are looking for small differences between large quantities. The solution to this problem 
is to modify stochastic perturbation theory so that it works directly with a boosted coupling. 

In Fig.|| we show the boosted perturbation series (from a 12 4 lattice) and compare it with the 
conventional series. Clearly the boosted series falls off much faster. Even though we are interested 
in regions where the boosted coupling is 50% to 100% larger the change is worth while. 

3. Perturbation Theory and Data 

We are now ready to compare perturbation theory with Monte Carlo data, to see if there is any 
remainder which we could attribute to a condensate. If there is, how does the condensate scale, like 
a 4 or a 2 ? 

In Fig.|| we show the result of summing the ordinary and the boosted perturbation series to 
various orders. We can see that boosted perturbation theory has a much better behaviour, it has 
almost 'converged' before we get to 12 loops, while the series in the original coupling is still 
changing visibly. Note that neither series shows any signs of diverging, though it is possible that if 
we go much further in loop number we will at some point see terms growing. 
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Figure 3: Convergence of the two expansions of the plaquette at j3 = 6.20. The series in the bare lattice 
coupling is shown with blue circles, the series in the boosted coupling with black squares. It can be seen that 
the boosted series has much better apparent convergence. 



The quantity which interests us here is the difference between perturbation theory and lattice 
Monte Carlo data, AP = P pert — Pmc- The horizontal red line in Fig.|| shows the value of the 
plaquette as measured in lattice simulations. Looking at the boosted perturbation series it does 
seem that P pert is measurably larger than Pmc- It would not be as easy to estimate the gap using the 
unimproved series, because the result would depend heavily on extrapolated values of the p n . 

To see whether AP scales like a 2 or a" we plot AP /a 1 against a 1 in Fig.|[ Most of the plaquette 
values are from [11], with some additional data calculated by QCDSF. If the non-perturbative 
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Figure 4: AP/a 2 plotted against a 2 , with both quantities measured in physical units. If AP scales like a 2 the 
data will be constant, if AP scales like a 4 the data will lie on a straight line passing through the origin. 
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contribution scales like a 2 , the AP/a 2 would be a constant, but if the non-perturbative contribution 
scales like a 4 then the data should lie on a straight line passing through the origin. We see that the 
data looks more like a dimension four condensate. 

We can use this graph to set a bound on any dimension-two condensate. Around a 2 ~ 0.05GeV~ 2 
we see that \AP/a 2 1 < 0.0015GeV 2 = (40MeV) 2 . Since the natural scale for a dimension-two con- 
densate is Aq CD , this is a rather tight bound. 

Using eq. ( |L2| ) we can find the dimension-four condensate from the slope of this graph, 

{^GG^i « 0.04(1) GeV 4 = (450 MeV) 4 . (3.1) 

This is about the right order of magnitude, though somewhat larger than phenomenological esti- 
mates [§, ||] which tend to lie in the range 0.012 GeV 4 to 0.024 GeV 4 . 



4. Conclusions 

Stochastic Perturbation theory gives long series for some simple quantities, allows us to test 
ideas that we could not look at with conventional perturbation theory. We see, for example, that 
boosted perturbation theory does dramatically accelerate series convergence. Although there are 
good reasons for expecting that perturbation theory diverges, we see no sign of this yet. Both the 
conventional and boosted perturbation series still look well-behaved. 

The gluon condensate, found from the difference P pert — Pmc, seems to be dimension four, not 
dimension two. The value of the gluon condensate is approximately 0.04(1) GeV 4 , but with large 
uncertainties. 
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